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If dark matter interacts with quarks or gluons, the mediator of these interactions is either directly 
accessible at the LHC or is so heavy that its effects are encoded in contact operators. We find that 
the self-consistency of a contact operator description at the LHC implies bounds on the mediator 
scale stronger than those found from missing energy searches. This translates into spin-independent 
elastic scattering cross-sections at a level < 10 -41 cm 2 , with direct implications for the DAMA, 
CoGeNT, and CRESST-II anomalies. We then carefully explore the potential of monojet searches 
in the light mediator limit, focusing on a Z' model with arbitrary couplings to quarks and dark 
matter. We find that the Tevatron data currently provides the most stringent bounds for dark 
matter and Z' masses below 100 GeV, and that these searches can constrain models for the DAMA, 
CoGeNT, and CRESST-II anomalies only if the mediator can decay to a pair of dark matter particles. 



I. INTRODUCTION 

If dark matter (DM) couples to quarks or gluons, its 
interactions can be discovered or constrained by LHC 
and Tevatron data in missing transverse energy (MET) 
searches. This has an interesting connection to direct de- 
tection experiments since the same couplings are probed. 
There is therefore the possibility that a one-to-one cor- 
respondence can be established between MET rates at 
colliders and scattering rates at direct detection exper- 
iments. Indeed, the status of DM direct detection, 
with the three experiments DAMA [I], CoGeNT 0, 
and CRESST-II [3] reporting excesses consistent with 
a 10 GeV DM candidate, cries out for an independent, 
astrophysics-free verification. This has recently stimu- 
lated significant work aiming at a consistent explanation 
of these experiments [H-iT?]. 

A rather general way of addressing this problem is the 
construction of an "effective theory" of DM that can be 
probed at colliders [T%rl2"7] . If the mediators of quark- 
DM or gluon-DM interactions are sufficiently heavy, then 
their effects can be safely integrated out and described by 
higher-dimensional contact operators. In direct detection 
experiments the elastic scattering of WIMPs on nuclei is 
well-described by a contact interaction so long as the me- 
diator mass is above a few MeV. In contrast, colliders like 
the LHC probe much higher momentum scales, and the 
validity of the effective description comes into question 
when the new physics is lighter than a few TeV. 

Here we will see that the very assumption that the DM 
interaction be consistently described by a contact oper- 
ator up to LHC energies implies perturbative unitarity 
bounds on the mediator scale that turn out to be much 
stronger than those imposed by MET constraints. Im- 
portantly, this assumption is also shown to be at odds 
with both the DAMA and the uncontaminated (surface 
event-free) CoGeNT regions. 
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The implication of this result is that if DAMA, Co- 
GeNT, and CRESST-II are indeed detecting DM, then 
the physics responsible for it must be directly accessible 
at the LHC and thus not generally amenable to a contact 
description. This gives hope to the prospect of confirm- 
ing or ruling out the direct detection anomalies in the 
near term. To this end we analyze in detail constraints 
arising from monojet searches at the Tevatron and the 
LHC in models where the mediator is a new massive vec- 
tor boson Z' with arbitrary couplings to the standard 
model quarks and the DM. 



II. UNITARITY BOUNDS ON CONTACT 
INTERACTIONS 

To quantitatively assess the range of validity of a con- 
tact interaction at the LHC, consider adding to the SM 
Lagrangian the operator 

o (i) 

where q is a SM quark and X the dark matter, assumed to 
be a Dirac fermion for definiteness. Clearly, this is a very 
simplified picture: a complete, unambiguous description 
of the theory necessarily contains additional operators. 
These are expected to become relevant before the energy 
probed E is such that the coupling E 2 / 'A 2 becomes large. 
By naive dimensional analysis, our simplified effective de- 
scription obtained by considering only the SM plus O is 
valid when A > E/in. 

Although this simple consideration already provides a 
useful determination of the range of validity of our con- 
tact interaction assumption, it relies on a somewhat ar- 
bitrary definition of "non-perturbativity." 

We can be slightly more precise regarding the valid- 
ity of the E/A expansion by requiring that our effec- 
tive operator description be consistent with perturbative 
unitarity. (Analogous arguments are familiar from pion 
scattering, and have similarly been used to find, e.g., 
an upper bound on the Higgs mass [28], the scale of 
fermion mass generation |29j . and on the dark matter 
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mass [30 .) In the presence of O, DM can be produced 
at the LHC via pp —> XX + where the dots stand 
for SM particles. The hard process contributing to these 
events is qq — > XX, and will be probed up to an energy 
E < 7 TeV. The amplitude for the process qq — > XX 
increases as a power of the center of mass energy y/s of 
the quark-antiquark system, and its probability eventu- 
ally exceeds unity, signaling a loss of confidence in the 
contact interaction assumption. The strongest energy 
dependence arises in the scattering of opposite helicity, 
color-singlet states. Specifically, the properly normalized 
initial and final states we consider are defined respec- 
tively as J2 a (\qfqJ) + \q%^)) / y/2N c ,jnth a = l,...,N c 
the color index, and (\XlXr) + \XrXl)) / 'y^, for which 
the (parton-level) amplitude reads 



M 



(2) 



with /3(s) = yl — /s the dark matter velocity. Up 
to a normalization, this amplitude coincides with the ze- 
roth order partial wave 
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Unitarity of the S'-matrix imposes the bound |ao(s)| < 
l/2/3(s) for elastic as well as inelastic processes (as a 
reference see, for example, [3T]). We thus see that O 
provides a sensible, effective description of DM-quark in- 
teractions only if 



A > 0.4 E/3 



2.(i TV [ -^L) 
7 TeV J 



(4) 



where /3 = /3(E 2 ) and E is the maximum energy flowing 
in the process. 

The physical interpretation of the above unitarity 
bound is that a UV-completion for O must appear be- 
fore Eq. Q is violated. This can manifest itself as new 
dynamics, or as higher order terms in the perturbative 
series. In either case our simplified description based only 
on O becomes unreliable when Eq. Q is violated, and 
a new, consistent theory of DM-quark interactions must 
be adopted. 

In calculations of pp scattering at the 7 TeV LHC one 
implicitly assumes a consistent and controllable descrip- 
tion of quark interactions up to ~ 7 TeV. 1 We therefore 
conservatively take E — 7 TeV in what follows. 

One may at first see this as an overly strong assump- 
tion since parton distribution functions (PDFs) will in- 
evitably suppress the characteristic scale probed by O 



1 In most searches, the kinematic cuts are typically of the order of 
a few hundred GeV, such that E ~ 7 TeV is justified. 
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FIG. 1: Here we compare the monojet rate coming from the 
contact operator O (cr contact) to that of a Z' model (cr) with 
the Z' coupling equally to quarks and dark matter. In both 
models we turn on only couplings to left-handed up quarks 
and dark matter and set rax = 10 GeV. In comparing the 
two models we fix A = m Z i /g, such that their low-energy 
predictions agree. The shaded, strongly coupled region g 2 > 
(3) 47r/v / 3 represents the unitarity constraint of the UV model. 



compared to the maximum center of mass energy of the 
LHC. However, although PDFs suppress the probability 
that highly energetic qq — > XX events are produced, 
they do not prohibit them. Thus the conclusion that new, 
UV-physics beyond O must appear if A < 2.6 TeV is not 
altered by the effect of PDFs. Instead, the relevance of 
PDF suppression becomes now a practical question: Is 
this new physics actually observable? How well can DM- 
quark processes be approximated at the 7 TeV LHC by 
a contact operator description that violates Eq. Q? 

These questions cannot be addressed unambiguously, 
since for example they inevitably depend on the specific 
UV-completion of the DM-quark interaction, the process 
under consideration, and the uncertainties inherent in 
the experiment. However, it is instructive to investigate 
under which conditions a tractable UV-completion for 
DM-quark interactions can efficiently be described by O 
in a given collider search. We focus on monojet plus MET 
events, and without loss of generality consider models 
where O is completed by the exchange of a resonance 
of mass gA. Now the unitarity bound is equivalent to a 
perturbative bound on the coupling, g 2 < 47r/v3, and 
the assumption A < 2.6 TeV is the statement that the 
mediator is kinematically accessible at the LHC. 

Clearly, if the mediator is exchanged at loop-level the 
coupling will be much smaller than unity and the reso- 
nance so light that a contact description is certainly in- 
adequate. Yet, we will now argue that even for tree-level 
completions the contact operator description is typically 
inappropriate when the unitarity bound is violated. 

There are in general two classes of tree-level comple- 
tions. The first consists of "s-channel" completions with 
interactions with a mediator <f> of the form qq/fi and XX<f>, 



3 




0.1 0.5 1.0 5.0 10.0 50.0 100.0 

m x (GeV) 



FIG. 2: Here we compare the direct detection upper bounds 
obtained from ATLAS monojets [231 [Ml EH (solid blue) and 
the unitarity constraint Q with E = 7 TeV (solid red line). 
We have also included the optimistic 5a LHC reach at 14 
TeV [23] (dashed blue) along with the unitarity bound, at 
E — 14 TeV (dashed red). All bounds are derived under the 
assumption that the quark-DM interaction O remains contact 
at LHC energies and with universal quark coupling. For ref- 
erence we include the DAMA 3a [2 [32], CoGeNT 90% CL g], 
and CRESST-II la and 2a preferred regions [3] as well as the 
90% CL XENON-100 bound [33]. 



and will be studied in detail in Section Mil The second 
class of models are "^-channel" completions in which the 
mediator is a color-triplet and has interactions of the 
form q<j)X. In both cases, as soon as Eq. ^ is violated 
the mediator becomes accessible at the LHC, and MET 
processes will tend to be dominated by resonance pro- 
duction (either qq — > <fig — > XXg or qg — > X<f> — > qXX 
depending on the particular completion), and the signal 
will typically be larger than that obtained from a contact 
interaction description. The resonance enhancement can 
be suppressed by pushing the theory towards the strong 
coupling limit g 2 ~ 47r/v3, but in so doing the theory 
becomes increasingly unpredictive. 

A quantitative measure of the impact of the resonance 
enhancement is shown in Fig. [TJ where we compare the 
monojet cross section of a contact description and an s- 
channel UV completion with the same value of A. From 
this we see that A above a few TeV is needed for an 
s-channel completion to produce a monojet signal com- 
patible with that of a contact interaction. Similar results 
can be found in [231155]. 

We therefore conclude that when A is smaller than a 
few TeV the observable predictions of a generic, tractable 
UV-completion and a contact operator will typically dif- 
fer greatly. A model-independent and conservative esti- 



mate of the bound is A < 2.6 TeV. 2 

The consistency constraint Eq. Q is much more strin- 
gent than the bounds obtained so far from missing energy 
searches [l3l[5l][53l[54l[26]. 3 This continues to be true at 
14 TeV, where unitarity requires A > 5.2 TeV, which is 
more stringent than the optimistic monojet bound with 
100 ftr 1 of data at the 14 TeV LHC [5J [55]. 

While in the above we chose to focus on the operator O, 
it is easy to see that one would obtain similar results for 
operators with different Lorentz structure, as well as for 
dark matter with different spin. Analogous constraints 
also apply to dark matter interactions with gluons. 

The consistency bound found above has important im- 
plications for direct searches of DM. In Fig.[2]we find that 
both the DAMA and uncontaminated CoGeNT regions 
are inconsistent with a contact operator description at 
the 7 TeV LHC, though the inclusion of surface events 
for CoGeNT may alter this conclusion [13 [33]. The di- 
rect detection implications of unitarity are even stronger 
in the case of velocity-dependent elastic scattering. In 
that case the bounds on the DM-nucleon scattering cross 
section in Fig. [2] scale down parametrically by powers of 
the DM halo velocity (v 2 ~ 10 -5 ). 

Thus, if DM is indeed the source of the DAMA and Co- 
GeNT anomalies then the mediator responsible for DM- 
quark interactions is already kinematically accessible at 
the LHC. 



III. DARK MATTER AND LIGHT Z'S AT 
HADRON COLLIDERS 

Once the contact interaction hypothesis is abandoned, 
matters become inevitably more model-dependent. Yet, 
irrespective of the details of the model, missing energy 
signals remain a characteristic signature of dark matter 
production at colliders. 

Jets and/or photons plus MET signals have been used 
before to bound dark matter [501 [5TJ [53J [53] and neutrino 
interactions [26]. Most of these efforts assume that the 
dark matter interacts with the SM via contact interac- 
tions. A first qualitative look at the effect of relaxing this 
assumption was presented in [5T] (see also [551 EH] ) i while 
a detailed analysis of the Tevatron monojet bounds for a 
specific model can be found in [55] , 

It would be useful to have a systematic and compre- 
hensive study of the light mediator limit that can easily 



2 In a previous version of the paper we attempted a different and 
more conservative estimate. We now believe that any such esti- 
mate would necessarily have limited applicability and we there- 
fore decided to take a more model- and process-independent ap- 
proach. 

3 Monojets currently provide a direct constraint on contact DM- 
quark interactions of order A > 500 GeV [20l ETtl l23l l2il [26] . 
Consistently with Fig. ^ such low values of A cannot be ac- 
counted for by tractable UV-completions; see also Fig. 7 of [24| 
and Fig. 7 of [26]. 
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FIG. 3: Upper bounds on y/glgx as a function of the DM 
mass when m Z i < 2mx, with m Z i <C 2mx (solid lines), 
and m z > — 100, 500 GeV (dashed lines). CDF provides the 
strongest bounds when the DM mass is less than about 30 
GeV. Note that the value of the mediator mass only matters 
near threshold mz' ~ 2mx, where the effect of the Z' width 
is also relevant. 



be applied to one's favorite model. Here we aim to make 
a step in this direction by analyzing in detail models in 
which the interaction between dark matter X and quarks 
q is mediated by a new vector boson Z' . 

We find that for light Z's and light DM the most con- 
straining MET search currently arises from the monojet 
data of CDF and the LHC: 

pp/pp -> j + MET. (5) 

In particular, we find the above monojet process to be 
more constraining than the analogous multijet + MET 
analyses. Moreover, the Tevatron's CDF data provide 
more stringent bounds than LHC data when the invariant 
mass of the DM system is small [26 , as shown in Fig. [3] 
These results are simply a consequence of the fact that 
CDF's monojet selection criteria are softer and thus more 
sensitive to a MET signal coming from low mass physics. 

A. Light Z's and MET signatures 

There are at least two well motivated ways to couple a 
new light Z' to the SM without affecting the huge phe- 
nomenological success of the SM [3"6"lf4"5] . The first is to 
couple the new boson via kinetic mixing with the hyper- 
charge vector. The second is to promote the baryon num- 
ber to a local symmetry. In either case the physics rele- 
vant here is effectively described by the SM Lagrangian 



plus £kin + Z' J^, where the first term contains pieces 
quadratic in Z' and X, while the second describes the 
interactions between the two sectors. In general 

^ = E^l<+.9 9 V)<Z + ^, (6) 

9 

where the precise structure of J x depends on the dark 
matter spin. Our MET bounds will be presented for 
fermionic dark matter, in which case 

J x = X^(g x+9x7 5 )X. (7) 

We will discuss how the bounds change for dark matter 
of different spin at the end of this section. 

(Note that even if the dark matter is assumed to be 
a SM singlet, loop effects can lead to couplings of X to 
the SM vector bosons of the order ~ 10~ 2 g. While these 
effects provide interesting model-dependent constraints 
on this scenario, they are irrelevant here.) 

MET rates at hadron colliders are set by the following 
differential cross section 

dcruET = 2J da q 9q9x > ( 8 ) 

where the smallness of the SM quark masses allows us to 
define the effective quark coupling g 2 = Iff^l 2 + |<?,f | 2 - An 
analogous expression for the effective dark matter charge 
is also a good approximation when mx < 500 GeV. In 
this range we can thus write g\ s» |#x| 2 + \Sx^- 

For each quark the full parameter space is therefore 
described by four quantities: mx , , s/9q9x ; an d the 
Z' width V z' ■ There are however two distinct physical 
regimes of interest where the number of parameters re- 
duces: light dark matter, 2mx < mz>; and heavy dark 
matter, 2mx > mz>. In the case of light dark matter, the 
mediator Z' can be produced on-shell, and subsequently 
decay to a pair of DM particles. Thus for a given mz 1 
MET searches provide an upper limit on the combination 

g q x yjBR(Z'^XX), (9) 

with a negligible dependence on the DM mass. (Here 
we have made use of the narrow- width approximation.) 
In contrast, for heavy dark matter the Z' intermediate 
state is necessarily off-shell and for a given dark matter 
mass, MET searches then probe ^Jg q gx , with a negligible 
dependence on mz> and Tz> ■ The intermediate regime 
2mx ~ mz> is sensitive to all four parameters. This 
is illustrated in Fig. [3] where the Z' mass dependence is 
apparent even in the heavy dark matter regime. In order 
to present our bounds in simple, two-dimensional plots 
we ignore this intermediate, somewhat fine-tuned regime 
in what follows. 

In the remainder we specialize to MET searches with 
a single jet in the final state, see ([5]), since we find these 
to be the most constraining. 
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FIG. 4: 95% CL monojet bounds on the Z' model from CDF and ATLAS in the low DM mass and heavy DM regimes (see the 
text for more details). The bounds are derived assuming that a single quark coupling at a time is switched on. In Eq. ( |10[ ) we 
show how to modify the bounds in the most general case. As a reference we also quote the bound arising from the assumption 
of universal coupling. The "plus" signs indicate the region were the CDF and ATLAS bounds cross each other (see also Fig. [3| . 



B. Monojet Bounds 

We simulate the signal using Madgraph_v5 46J, then 
apply the selection criteria of CDF |47rf5D] and AT- 
LAS [5JJ using Pythia 6 . 4 [55] and cluster with Fast j et 
2.4.4 [53] . 

In Fig. [i] we present the 95% CL bounds on DM-quark 
interactions in both the light (left) and heavy (right) DM 
regimes. The limits are derived assuming that only one 
quark coupling at a time is switched on. When several 
couplings g q are turned on our bounds read 



E 



9q 



.abound 

y q 



< 1, 



(10) 



with n = 2(4) for the light (heavy) DM regime and g q ound 
being the bound on the quark flavor q shown in Fig. [4] 
The resulting bounds under the universal quark coupling 
assumption are shown as dashed black curves. Notice 
that the up and down quark bounds are stronger than 
the curves corresponding to the sea quarks as a result of 
PDF enhancement. The curves asymptote to a constant 
when the DM invariant mass is small compared to the 
Pt cut on the leading jet. 

At large DM invariant mass, the veryHighPT cut from 
ATLAS [ST] sets the strongest limit, whereas for low DM 
invariant mass the CDF 1 fb _1 analysis [49] wins. We 
note that although a CDF 6.7 fb _1 analysis based on a 
jet Et analysis exists [50) . it is less constraining than the 
earlier 1 fb _1 counting analysis. In Fig. [4] we have taken 



care to plot only the strongest of the collider monojet 
bounds. The point at which the crossover occurs is shown 
as a "plus" in the figure. Note that the mass at which 
the Tevatron starts setting the dominant bound is larger 
for quark flavors with a stronger PDF enhancement. 

Although we have focused on fermionic DM thus far, 
it is straightforward to translate the bounds to DM with 
other spins. First, note that the left plot of Fig. [4] applies 
to any DM candidate so long as Z' —> XX is allowed. For 
light DM, the bounds are roughly rescaled by the spin 
degree of freedom. For example, for a complex scalar 
DM particle the bounds on y/gxgq in the right plot of 
Fig. [4] are decreased by a factor l/\/2. (Here gx stands 
for the coupling of the XX Z' cubic vertex.) 

Lastly, we can translate our monojet bounds to direct 
detection cross section limits. As is visible from Fig. [4] 
the constraints on heavy DM are quite weak. For ex- 
ample, in the case of universally coupled Z' with a 10 
GeV DM mass, the spin- and velocity-independent DM- 
nucleon elastic scattering cross section is bounded by 



&nX 



< 10 



" 34 cm 2 



20 GeV V 
m z > ) 



(11) 



Note that in this regime the upper bound on the Z' 
mass, mz> < 20 GeV, implies that the constraint be- 
comes weaker as the Z' mass decreases. In contrast, when 
mz> > 2nix = 20 GeV the limits are stronger, especially 
for narrow Z's, implying 



OnX 



< 10~ 36 cm 2 



20 GeV 



Tz'/mz' 

10- 2 



(12) 
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Now the bound for fixed Tz> / mz 1 becomes more strin- 
gent as the Z' mass increases, approaching a n x ^5 10~ 40 
cm 2 when mz' ~ 500 GeV. 

IV. CONCLUSIONS 

We have seen that the very assumption that the DM- 
quark interaction remains contact at the LHC provides 
strong constraints on the scale of new physics. Specif- 
ically, perturbative unitarity arguments force the scale 
suppressing the higher dimensional operator to be above 
a few TeV. These bounds are notably stronger than those 
obtained from monojet searches at the LHC. 

In light of these results, the DAM A favored region and 
the (uncontaminated) CoGeNT 90% CL signal region 
cannot be explained by interactions that remain contact 
at the 7 TeV LHC. 

We therefore relaxed the contact interaction assump- 
tion and analyzed in detail the constraints arising from 
monojet searches at the Tevatron and the LHC. We fo- 
cused on models in which the DM-quark interaction is 
mediated by a light Z' with generic couplings to the SM 
quarks and the DM, and presented the bounds in a way 
that our results are readily applicable to the reader's fa- 
vorite scenario. 

We find that the strongest bounds are currently pro- 
vided by the Tevatron's CDF experiment for DM and Z' 
masses below ~ 100 GeV, with the LHC taking over for 
heavier masses. 

Translating our monojet bounds to limits on the DM- 
nucleon elastic scattering cross section we find that these 
searches can constrain models for the DAMA, CoGeNT, 
and CRESST-II excesses only if Z' — > XX is allowed and 
if the Z' is a narrow resonance. 



Note Added 

After the submission of our arXiv paper, the preprint 
|54) appeared with a discussion claiming that the unitar- 
ity constraint of the present paper is in practice unim- 
portant for rather low cutoffs. This conclusion is based 
on the inspection of the XX invariant mass distribution 
for a MET search, in the case of a contact operator with 
A ~ 600 GeV. From this the authors find that the num- 
ber of events that "violate unitarity" (i.e., those events 
with m,x X > A/0.4 w 1.5 TeV) is negligible for light 
dark matter. However, as we clarified above, the unitar- 
ity bound implies that a consistent description requires 
that new physics appear at invariant masses below 1.5 
TeV. A more telling measure of the validity of the con- 
tact operator description would be a comparison between 
a contact operator and a UV-complete theory leading to 
the same A. This comparison was already made in [53] 
for monojets, where it was shown that in order for an s- 
channel completion to agree with the contact interaction 
prediction one should have couplings of order g ~ 5 — 10, 
in agreement with our Fig. [T] With such large couplings, 
yet higher dimensional operators are expected to become 
relevant and the monojet prediction is plagued by 0(1) 
uncertainties. We interpret this as a confirmation of our 
conclusions. 
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